Drastic effects of damping mechanisms on the third-order optical nonlinearity 
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We have investigated the optical response of superradiant atoms, which undergoes three different 
damping mechanisms: radiative dissipation (7r), dephasing (7d), and nonradiative dissipation (7n). 
Whereas the roles of 7d and -y n are equivalent in the linear susceptibility x > the third-order 
nonlinear susceptibility \^ drastically depends on the ratio of 7d and ■y n : When 7d <C 7n, x' 3 ' 
is essentially that of a single atom. Contrarily, in the opposite case of 7d 3> 7n, x' 3 ' suffers the 
size-enhancement effect and becomes proportional to the system size. 

^ [ PACS numbers: 42.65.-k, 42.50.Fx 

o : 

There has been much interest in optical responses of finite-sized systems. When the eigenstates of the system arc 
delocalized in space, the transition dipolc moment between the ground state and the lowest excited state scales as 
iV 1 / 2 , where N is a parameter representing the system size Q. This N 1 / 2 scaling rule is the origin of unique optical 
\ responses of finite-sized systems, such as size-enhancement of third-order optical susceptibility, x^ 0, IE El IE IE 01 • 
In conventional theories on the size-effects of optical responses, whereas the size-dependences on eigenenergies and 
transition dipole moments were carefully taken into account, damping effects were often treated rather crudely by 
simply introducing phcnomcnological damping constants. However, because the system size would affect not only 
t-H ' eigenenergies and transition dipolc moments but also damping rates, more rigorous treatment on the damping effects 
are desired. Particularly, when the material is irradiated by resonant light fields, the magnitude of the optical 
nonlinearity is strongly sensitive to the damping rates. It is therefore indispensable to exclude phenomenology on the 
damping effects for quantitative evaluation of optical nonlinearity. 

Here, targeting solid-state nonlinear optical devises in mind, we investigate the optical response of finite sized sys- 
I | terns under a situation where the system suffers three different damping mechanisms: radiative dissipation, dephasing, 
. and nonradiative dissipation. The latter two damping mechanisms are brought about by coupling to environmen- 
tal degrees of freedom such as phonons. The effects of nonradiative dampings have been considered in detail in 
Ref. |fj, but the damping constants are introduced by hand, independently of the model. (The radiative damping is 
incorporated through the self-consistent Maxwell fields in their formalism 0.) Spano et. al. pioneered the theories 
without phenomenology on the damping effects, where the damping dynamics of the system is explicitly defined in 
the model 3,0. As for nonradiative damping effects, they introduced the homogeneous dephasing alone and no dis- 
sipation was explicitly treated. Although particular aspects of damping effects on the size-dependence of nonlinearity 
have been considered, the interplay of different damping mechanisms is still a subject of importance. To reveal the 
overall effects brought about by three different damping mechanisms, it is indispensable to treat the radiative and 
nonradiative dampings on equal footing and to make a clear distinction between dephasing and nonradiative dissipa- 
tion. It is shown here that the third-order nonlinear response is drastically dependent on the ratio of two kinds of 
nonradiative damping rates: the ratio determines whether the size-enhancement of nonlinear response occurs or not. 



OO 



■ This fact implies that particular damping conditions could provide a novel resource for size-enhancement of nonlinear 
response. 

The objective of this study is to investigate the third-order nonlinear optical response, taking account of different 
damping mechanisms explicitly. As a simplest model of a nonlinear optical system with finite size, we consider a 
superradiant system composed by N identical two-level systems (hereafter referred to as "atoms") with transition 
frequency f2 [8j, which suffers, individually at each atom, both dephasing and nonradiative dissipation. Such damping 
mechanisms become particularly significant if the atoms are embedded in a solid-state environment, e. g., quantum 
dots in a microcavity. The equation of motion for the density matrix p of atoms is the superradiant master equation 0, 
IE IE] supplemented with the terms describing dephasing and nonradiative dissipation |10| . It is given, omitting h and 
\x (transition dipole moment) for notational simplicity, by 

% = -i[n + n- mt (t), P \ + (c I + c d + c n ) P , (i) 

at 
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Ho = Y,to*W ( 2 ) 
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Wta*(t) - E(t)S* + E*(t)S, (3) 
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£ r p = -|(2Sp£t - S^Sp - pS^S), (4) 



7d 
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where and Sj are the Pauli creation and annihilation operators at jth atom, and S — y~] ■ Sj is the collective operator. 
-E'(i) represents the positive frequency part of the applied electric field, and the rotating wave approximation is used 
in Ti.i n t(t). 7 r , 7d and 7 n represent the single-atom rates of radiative decay, dephasing, and nonradiative dissipation, 
respectively. It is of note that the atoms interact with the electromagnetic field via the collective operator [see Eqs. © 
and whereas dephasing and nonradiative dissipations occurs independently in each atom [see Eqs. (J5J and J^J]. 

Based on this model, we investigate the linear and the third-order nonlinear optical responses. To the end of 
investigating up to third-order response, we are concerned with the following expectation values: (sj), (s{Sj), (s^Sj), 
and (s]sjSk), where the expectation value of an operator A is given by (A) = Tr{pA}. Remembering the fact that all 
atoms are equivalent, the number of independent variables are greatly reduced. For example, it is apparent that (sj) 
is independent of the site index i. We use the following notations: 

(*> = (*>> (7) 
(siSj) = I ,° ^ rl J l , (8) 

(sU) - { H*) ( * =j) f9l 
[ 4 3> ~ \ ( S t s '} (i ^ j) ' W 

f (j = k) 

(slsjSk) = < (sUs) (j ± k, i = j or k) . (10) 
l(sW> (i^j,j^k,k^i) 

The equation of motion for (A) is given, using Eq. (JTJ, by ^(A) = -i{[A,H + H in t{t)}) + (j I /2)({S +: A]S- + 

s+[a, S-]) - ( 7d /2) EM A ' a j*jl a i s A) + (^/ 2 ) EM> A fa + 4^> Dennin S T -Ac by 

r o ,6,c = (07r + Hi + C7n)/2, (11) 

the equations of motion for (s) etc are given as follows: 



d( Sl )/dt = {-m-r NAA )( Sl )-iE, (12) 

d(ss)/dt = (-2in-r 2J v_a,a,2)(*a> -2iE(s 1 >, (13) 

d(sU)/dt = (i£*(s 1 ) + c.c.)-r 2)0 , 2 < s t s )-r 2JV _ 2 , 0! o(s t s / ), (14) 

d^s')/dt = (i£?*<ai)+cc.)-r 2i o,o<a t a)-r 2JV _ 2i2 , 2 (at fl '>, (15) 

d<a 3 )M = (-in-r w ,i 1 i)(s3>+2LB(ata)+r 2W _ 2j o,o<a t aa), (16) 

d{sUs)/dt = (-in- r JV+2 ,i,3){s t ss) -r 2A r_ 4 ,o,o(s t ss') + is*( ss )-i£;((/ s ) + ( s t s ')) j (17 ) 

d{sUs')/dt = (-\tt-T aN - 6 ^){sUs') -T Afl . {shs) +iE*{ss) -2iE(sU'), (18) 



where (si) and (S3) denotes the first- and third-order components of (s). Although not explicitly indicated, (ss), (s^s) 
and (s^s') ((s^ss) and (s^ss')) in the above equations are the second- (third-)order quantities. It is of note that, in 
the above equations of motion, the dependence on the system size N appears only through the enhancement of 7 r . 

We can easily obtain the stationary solutions of these simultaneous equations. Assuming that E(t) is monochromatic 
as E(t) ~ e~ iuJt , and introducing / ,6 jC (w) by 

/«,6,c(w) = (w-n + ir , 6)C )- 1 , (19) 

(si), (ss), (s^s) and (s^s') are given as follows: 
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(si) = f N ,x,i(u>)E, (20) 
(ss) = f N ,x,i{u)f N -x^ x {u)E 2 , (21) 

\ _ iV7r + 7d + 7n ( 7d + 7n A I f / x|2|p|2 / 99 ^ 

Eq. (1221 shows that, in determining time-independent quantities such as (s^s) and (s^s 1 ), the ratio of damping 
constants play crucial roles. This feature was also observed in conventional theories on the nonlinear susceptibilities 
using phenomenological damping constants 0, ld> l9- I n terms of the second-order quantities, (s^ss) and (s^ss 1 ) are 
given by 



(shs) \ = ( f N +2 i 3M ir 2JV-4,o,o \ / -E*(ss) + E(sKs) + E(sW) 
(sUs')J { ir 4 , , / 3 l 633 Mj ^ -E*(ss) + 2E(sW) 



(23) 



and (53) is given, in terms of (s's) and (s^ss), by 

(«s> = /jv,i,i(w)[-2 J B( S t s ) + ir 2A r_2 1 o,o(* t fl*>]. (24) 
Thus, we obtain the linear and third-order susceptibilities per one atom as follows: 

X (1) M = y = /«,i,iH (25) 



(sis) 

~2 , „ lo + ll2jV-2,0 : 07 



IE 1 ! 2 ZJV - Z ' U ^'|£;|2 S 



(26) 



both of which are free from phenomenological treatment on the damping effects. In the following part of this study, 
we discuss how x^ depends on the relaxation parameters (7,., 7d, 7 n ) and the number N of atoms. Recent nanotech- 
nologies aim to fabricate clean quantum systems with long coherence times, and extensive efforts on suppressing 7d 
and 7„ are being made. In the following part of this study, we restrict our attention to a case where 7d and 7 n are 
well suppressed and satisfy (7d,7n) <C 7r- 

Firstly, we discuss the limiting case of 7d — > 0. In this limit, it is easily confirmed that (s^s) = (s's'}(= 
\fN,o,i{w)\ 2 \E\ 2 ) and (s^ss) = (s' f ss')(= |/ A r i0 ,i(w)| 2 /jv-i,o,i( w )l^'| 2 £')- These equalities imply that N atoms respond 
to the electric field cooperatively, as a spin N/2 object, x is reduced to the following form: 

x£LoM = -2|/Ar, , 1 (a ) )| 2 / A r, , 1 ( W )/ Ar _ li0 , 1 ( W )/- 1 il ^). (27) 

(3) 

This equation reveals that xXJ-to depends on the system size N only through the enhancement of 7 r . In the off- 
resonant frequency regions, xl^— >o — — 2/(w — fi) 3 , which is independent of N. Thus, in the limit of 7d — > 0, the 
optical nonlinearity is essentially that of a single atom, except for minor corrections around the resonant frequency 
region. 

Next, we consider a more general case of 7d ^ 0. We should remark that, when (7d,7n) *C 7r is satisfied, (s^s) and 
(s's 1 ) are reduced to the following forms: 

, t , N~f d + Njn m2| F |2 / osA 

< s s > " 7d + JV7n \ e \ (28) 

( S V) ~ \f Nxl (")\ 2 \E\ 2 (29) 

7d + N'fn 

In case of 7d <C 7n, all of (ss), (s^s), and (s^s 1 ) are of the same order (~ |/| 2 |i?| 2 ), and the nonlinear susceptibility 
is given by Eq. 1(771) . Contrarily, in the opposite case of 7d 3> J n , (s^s) becomes much larger than (ss) and (s^s'). 
((sU) ~ iV|/| 2 |£| 2 , whereas (ss) and (sV) ~ |/| 2 |#| 2 .) Then, Eqs. © and © suggest that (shs), (sW) and 
X^ become almost proportional to (s^s). Using the fact that (s^s) is magnified by a factor (Njd + N~/ n )/(jd + N~/ n ) 
in comparison with the 7d — ► case, we obtain the following approximate expression of the third-order nonlinear 
susceptibility: 

~(3)/ \ Nj d +Nj n (3) , , ,„ n , 
X M " 7d + iV 7 n X ^° M - (30) 
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FIG. 1: Comparison of x^l and its approximate value |x' 3 '|, when N = 5. The solid, dotted and broken lines show x' 3 " 1 f° r 
{jd/ji, 7a/ji)=(0, 0.1), (0.05, 0.05) and (0.1, 0). The thin dotted lines show x' 3 ' with the same parameters. Note that x' 3 ' 
coincides with x for (7d/7r, 7n/7r)=(0, 0.1). 
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FIG. 2: Temporal behavior of Ix' 3 ^! assuming that irradiation of monochromatic electric field, E{t) = Ee~ lujt , starts at t = 0. 
The following parameters are used: N = 5, u> = f2 + 57 r , 7n = 0, and 7d/7r = (solid line), 0.01 (dotted line), and 0.05 (broken 
line). Transition to the enhanced value takes place at t ~ iV/(27d). 

In Fig.Q] the approximate susceptibility \^ is compared with the rigorous susceptibility x*- 3 -* ■ The figure demonstrates 
that x^ serves as a good approximation of ■ 

Now we discuss the implications of Eq. (|30|l . As far as the linear optical response is questioned, Eq. (|25|l indicates 
that the roles of dephasing (7^) and nonradiative dissipation (7„) are equivalent in determining the linear optical 
response. In contrast, when one questions the nonlinear optical response, the roles of two damping mechanisms are 
no more equivalent: The prefactor of the RHS of Eq. (|3U1) indicates that magnitude of x is sensitive to the ratio 
7d/7n, even if both 7,3 and 7 n are much smaller than j T . When 7d <C 7 n , X^ 3 "* is essentially independent of the system 
size N. Contrarily, when 7d 3> 7n, X suffers the size-enhancement effect. This observation demonstrates that it is 
indispensable for quantitative evaluation of nonlinear susceptibility to discriminate two damping mechanisms and to 
treat them non-phenomenologically. 

One might feel uneasy about the fact that x^ 3 "* is indefinite at 7d = 7n = 0. In order to resolve this problem, we 
investigate the transient optical response by considering a situation where monochromatic field, E(t) = Ee~ lu)t , is 
switched on at t — 0. By inspecting Eqs. (|14l) and (|15l) . we can find that x has two relaxation rates, Nj T (= t^ 1 ) 
and 27d/iV + 2j n (— t^ 1 ), before attaining to its stationary value. The temporal behavior of Ix^l is numerically 
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pursued in Fig- El The figure clarifies that x^ first relaxes to the unenhanced value within a short time (t < Ti), and 
the size-enhancement effect emerges gradually in the later stage (t ~ t 2 ). When 7d = 7n = (solid line in Fig. 
the size-enhancement does not take place, and x^ = X^— >o forever. More generally, when one is concerned with a 
transient behavior (t < T2), the size-enhancement effect is not expected; it is expected only after a long time, t > T2. 

Finally, we comment on the relevance to previous studies. Mathematically, x should be evaluated by the station- 
ary solutions of equations of motion for (sj), (s\sj) 1 etc @. Because these quantities generally have dependence on 
the site index i, it is usually difficult to obtain analytic expression of x^ ■ (By taking the eigenstates of the system 
Hamiltonian TYo as the basis, one may diagonalize the unitary part of Eq. Q). However, the basis generally does 
not diagonalize the damping part of Eq. £Q) simultaneously 4]. The conventional expansion for x 1 1 1 1] is obtained 
by approximately neglecting the off-diagonal part.) In the model of our study, an analytic form of x' 3 ^ is obtained 
without any approximation by virtue of symmetry of the system, and it is revealed that tiny difference in the damping 
rates in Eqs. (|14() and l|15fl may result in drastically different optical response, as observed in Fig. [21 A novel prediction 
in the present study is that the size-enhancement may take place, even when there is no transfer of excitations among 
the atoms. 

In summary, we have analyzed the third-order nonlinear susceptibility of superradiant atoms, which undergoes 
three different damping mechanisms: radiative dissipation (71), dephasing (7a), and nonradiative dissipation (7 n )- 
The analysis is based on the superradiant master equation supplemented with effects of dephasing and nonradiative 
dissipation [Eqs.JU-©]- The linear susceptibility x an d the third-order susceptibility x^ P er one atom are given 
by Eqs. I|25l) and l|26|). and x^ is well approximated by Eq. i|30|) . Whereas the roles of 7d and 7 n are equivalent 
in x^ [see Eq. (|25|) ]. they are no more equivalent in x^ [see Eq. (|30|) ]: x'" 3 '' depends on the ratio 7d/7n- When 
7d <C 7n, x'" 3 '' is essentially that of a single atom. Contrarily, when 7d ^ 7n, suffers the size-enhancement effect 
and becomes proportional to the system size N. These observations indicate that, for qualitative evaluation of x^ 3 \ 
it is indispensable to distinguish 7d and 7 n clearly, and to handle them in a non-phenomenological manner. 
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